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O ' Abstract 

o : 

We consider a linear partial integro-differential equation that 
arises in the modeling of various physical and biological processes. 
^ , We study the problem in a spatial periodic domain. We analyze 

numerical stability and numerical convergence of a one step approxi- 
mation of the problem with smooth and non-smooth initial functions. 

Keywords: stability; convergence; smoothness; full discrete problem; infi- 
nite domain; periodic domain. 

1 Introduction 

Study in convolution model of phase transitions (initial value problem) is 
of ongoing interest. Many scientific problems have been modeled by reac- 
tion diffusion equation and advection reaction diffusion equations. A lot of 
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models consider nonlocal type diffusion operators [21 [31 HI El HUl [13]. In this 
article we analyze stability, accuracy and rate of convergence of a simple 
approximation of such a linear partial integro-differential equation (IDE) 

ut{x, t) = e J{x - y) {u{y, t) - u{x, t)) dy = ILu, (1) 

Jn 

which is the linear part of the IDE [21 [31 [SI [3 [SI [IHl [IB [12] 



ut = e{ J{x-y)u{y,t)dy-u{x,t) J{x - y)dy \ + f{u). (2) 

Here the initial condition u{x,0) = Uq{x),x G Q where f2 C M, /(«) is a 
bistable nonlinearity for the associated ordinary differential equation Ut = 
f{u) and J{x — y) is a kernel that measures interaction between particles at 
positions x and at y. Here it is assumed that the effect of close neighbours 
X and y is greater than that from more distant ones; the spatial variation 
is incorporated in J. We assume that J is a non-negative function satisfies 
smoothness, symmetry and decay conditions. 

We present the stability and convergence of the rescaled integro-differential 
equation model in a periodic domain. Here we consider spacial [0, 1] peri- 
odic domain. For simplicity of notations, from here J, u mean functions 
in the periodic domain and J°°, u°° mean functions defined in the infinite 
domain. If we choose spatially one-periodic initial data u{x, 0), then for all 
a; G M and t G M+ 

u{x,t) = u{x + 1, t). 

Then a with kernel function J°°{x), ([T|) can be written as 
Ut = e J°°{x-y){u{y,t)-u{x,t))dy 

= e J{x — z) {u{z,t) — u{x,t)) dz, (3) 

Jo 



where 



J{x) = J2 J'^ix-r) (4) 



and X G [0, 1]. Two sample kernels are shown in Figure [H Here we consider 
e = 1 only since it does not affect our analysis. 

In [3] , the author study stability and accuracy of an one step approxima- 
tion the IDE considering infinite spatial domain. They first write a discrete 
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Figure 1: Examples of kernel functions in [0, 1] periodic domain defind by 
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equivalent of the integral equation, then use the forward Euler method for 
time discretisation. They present some numerical results to demonstrate 
the rate of convergence for smooth and non-smooth initial functions. 

In |10], Duncan et. al. consider ([2]) in a spatial periodic [0,1] and 
[0, 1]^ domains. They approximate the problem using piecewise constant 
basis function with collocation and mid point quadrature rule for space 
discretisation, then they use some standard ode solver for time integration. 
They present some numerical results to demonstrate their scheme. In [8], 
the author consider several local and nonlocal operators that contains ([T]). 
He approximate the models using finite difference schemes and discuss some 
stability issues. 

Here we consider the problem ([3]) in a spatial periodic domain, and 
study stability and accuracy of a simple scheme using the Fourier series and 
the discrete Fourier transform definitions. We Organize our article in the 
following way. We present a simple space time discretisation in Section 2 
followed by stability analysis of the scheme in Section 3. We analyze accu- 
racy of the scheme in Section 4 considering smooth initial function, whereas 
in Section 5, we discuss the same considering non-smooth initial function. 
We study accuracy of the semi-discrete time dependent scheme in Section 6. 



We finish our study in Section 7 with some numerical experiments and dis- 
cussions. 

2 Numerical approximation 

We approximate ([T]) as 

du,{t) ""-' 



dt 



hJ2j{xj-Xk){Uk-U,) (5) 

fe=0 

for each j = 0, 1, 2, ■ ■ ■ , A^ — 1 where Uj{t) ^ u{xj,t) and Xj = jh. Applying 
Euler's method we can again approximate the semi-discrete problem (^ as 

N-l 

fc=0 

(Af-1 \ Af-1 

l-hAtJ2 Ji^j -Xr)] + hAt J2 Ji^j - ^r)U^- (6) 
r=0 / r-=0 

We use Fourier series and discrete Fourier transforms throughout. Be- 
fore the main discussion let us introduce some necessary definitions and 
theorems. The Discrete Fourier Transform for a periodic function (DFT) 
can be defined as 

Af-1 
3=0 

where u{x,t) is a periodic function with period 1 and its inverse Fourier 
transform is defined as 



N-l f 



Uj 



hJ^Uke'^'"''"' or Uj = h J2 Ukc'^'"''"'- 

k=0 fc=-^+l 



For any function periodic function f{x) with period 1, its Fourier series can 
be defined as 

fix) = Yl f^"^^^^^ ^^^^^ ^- = / /(^)e"^'"""rfx. 

n=-oo J^ 

For any 1-periodic complex-valued function u{x) the relation 

oo „i 

, io 



where «„ are the Fourier coefficients of u{x) and for a real- valued function 



uix] 



^2 1 °° /"l 



u^{x)dx 



N-l 


Af-l 






E -/ = 


■ aE^'^ 


or 


"sL-,, 


j=0 


k=0 







is called Parseval's relation for the Fourier series where a„, &„ are Fourier 
coefficients. For the Discrete Fourier Transform of u(x), the relation 



\Uk\ 



is called Parseval's relation where Uk is the DFT of Uj, both have the same 
length A^. 

Now we will discuss the relation between the DFT and the Fourier coef- 
ficients. From the DFT definition Uk = h^-~Q u{jh)e~^'^'"^^\ the Fourier 
coefficients are 

Uk= f u{x)e~'''^^^dx 
Jo 



and the corresponding Fourier series is 



fc=— oo 



Thus 



u>= = '"E E 

j=0 r=—oo 



N-l oo 



M„e'""""-'e ■' 



oo 



= h^ UrN6N{r -k) = ^ Uk+mN (7) 

r=— oo m,=— oo 

where r — k = mN and 

J 1 j = mN for some m E 7j 
1 otherwise. 

The relation ([7]) is known as the discrete Poisson sum formula [15]. The 
Fourier coefficients of the kernel function in [0,1] are given by 



pi /'OO 

Jo J -oo 



(9) 



is the relation between the Fourier coefficient of J(x) defined in (jl]) and the 
continuous Fourier transform of J°°{x). 

Multiplying both sides of ([6]) by }ie~^'^'^^^^ and summing over j 

N-l N-l / N-1 

j=0 i=0 \ r=0 

Af-1 / N-l 

+h Y, e-^^-'^'^f/; hAt Y J{ 

r=0 \ i=0 






That is, 

r=0 j=0 

= gih,At,k)U^ 

and clearly 

f/," = ^"(/i,At,fc)t/°, (10) 

where 



g{h, At, k) = l + hAtY Ji^j - ^r) (e"*^"'^^^ - l) = 1 + At (jk - Jo) 

3=0 

since J is 1-periodic. 



3 Stability analysis 

To show the stability of the scheme we need some reasonable restrictions on 
J{x). We will examine properties of the DFT of J{x) under some reasonable 
hypothesis on the function J{x). We use the bounds obtained below to get 
a bound on g{h, At, k). 

Lemma 1. Assume that 

Al. J(x) > 0. 

A2. J{x) = J{-x). 

A3. /^ J{x)dx = 1. 



A4. Aj(a;)<o/ora:G(0,|) 

A5. ifc>0/or-f + 1< A;< f 

then < Jo — J/t < 2, Jo — Jfe < 2 where — Y + l<fc<Y. Furthermore, 
Jfc>0. 

Proof. We have 

J;. = /i ^ J(x^)e-*'=2"^'- and Jo = /i ^ J(x,). 

r=0 r=0 

So 

Jo-Jk = hY,Axr){l-e"'''^^^^). 

r=0 

For simplicity from here we take A^ to be even. Using the symmetry of J(x) 
in [0,1] 

N_ 
2 



Jo — Jfe = 2/i ^^ J(xr-) (1 — cos(A;27rXr)) > 

r=0 

as /i > and J > 0. Also 1 — cos{k2TcXr) < 2. Thus 

— 1 

Jq- jj^ = Ah^J{xr) < 4 / ' J(x)Jx = 2 / J(x)Jx = 2. 

So we conclude < Jo — Jfc < 2. The result Jq — Jk ^ '^ follows from the 
definition of J and by using similar steps to those of [1]. D 

Lemma 2. If J{x) satisfies Al - A3 then 1 > g{h,At,k) > 1 - CAt for 
some C > and if J{x) satisfies A4 - A5 from LemmaUl o,s well then 
C = 2 and 

\g(h,At,k)\<l for all < At < At* = ^. 

Proof. Proof of this Lemma follows directly from Lemma [1] D 

The stability result follows from the following theorem. 

Theorem 1. If J{x) is a periodic function in [0, 1] and satisfies Al - A5, 
then there exists At* > given by Lemma\^such that ||f/"||/i < ll^'^IU for 
all < At < At* and n > 0. 



Proof. We have 



n-l Af-1 2 ^^1 



j=0 j=0 j=0 



N-l N-l 

\tt0\2 _ iir; ,, 

j II /i) 



< h'Y,\U!\'<hJ2\Uj\' = \\U^'' 



j=0 j=0 

which gives the result. D 

4 Convergence analysis of the fully discrete 
approximation 

Let the Fourier series of u{x, t) be 

oo 

,t)= Y,u,{t)e^^'- (11) 

where 



u(x, 

j=-oo 



Jo 
and let the Fourier series expansion of J{x) be 



where 



Substituting ([II])-([l2D in (^ 



J{x) = 


>: 

j=—oo 


j,^ij2.x 


Jo 


J{x)e 


-'^^^^dx. 



(12) 



oo 



d 

d 



Y^ 7%(^)e'''"" = Y^ %(t) J,e*^'2-- - ^ n,(t)e'^'2"" 

j=—oo \j=—oo j=—oo 



^^Uj{t) = {Jj - Jo)Uj{t) = QjUjit) 

with Qj = {Jj — Jo) where Jj is the jth Fourier coefficient of the kernel 
function J{x). Solving the above equation we have 

u,it) = e**%(0), 



where 

Thus the exact solution of the IDE ([T]) can be written as 

oo 

«(x,t)= ^ %(0)e«^V^'2"^ (13) 

j=-oo 

Now taking the inverse of the discrete Fourier transform in (11 Op we get the 
approximate solution of ([T]) as 

JV 
2 

U^= Y^ g''{h,M,k)Ule'^^^^K (14) 

From ([13]) and ([HD 

JV 
oo 2 

u{xj,t^)-UJ' = ^ Mfc(0)e'^'=*'"e'2"'="^ - ^ ^'"(/i, At, A;)[/°e''='"^^ 

JV 
2 



J] (wfc(0)e«'=*™ - r (/i, At, A:)[/o) e' 



|fc|>f 
Now taking the inner product on flT^ and applying Parseval's relation 



JV 
2 



|«(x„U-t/f||^ < Y. \ume''^''--9'^{KM,k)Ul 



fc=-f+i 



+ Y |wfe(0)e^''*™| . (16) 

|fc|>f 

The Poisson summation formula gives 



'^k — /_^ '^k+rN — 7 , 



oo 

oo,m. 



'^k+rN 



where "u^^jy ^^ ^^^ coefficient of the Fourier series at tm for the periodic 
function u whereas u^^^q represents the CFT for the nonperiodic (infinite- 
dimensional) case. A calculation similar to that leading to ([9]) gives 

{i^(0) = / u{x, 0)e-*'=2""dx = ul = it^ {2Tik) = w^'°. (17) 

Jo 
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So the first part of tlie riglit-liand side of fll6p can be written as 



2 



J2 Mfe(0)e^^*'" - ^™[/^° < J2 Ke^^*'"-^™) 



Wt 



',01 2 



+ 



JV 
2 



fc=-T + l 



^'"(/i, At, A;) J2 ^°^'° (27r(fc + rN)) 



Before giving the main convergence results let us introduce a lemma with 
reasonable restrictions on J{x). 

Lemma 3. // J(x) is a nonnegative even periodic function in [0, 1] and 
monotone- decreasing in [0, |] then (jj < and \Jo — Jj\ < | where Qj = 
{Jj - Jo) with Jj = J^ J{x)e-'^^^^dx. 

Proof. The proof of Lemma [3] follows directly from the definition of jj. D 



Now 



|e*^*| < 1 and \g{h, At,k)\ < 1. 



So 



|e^^*'"-^'"(/i,At,A;)| = | (e^^^*)™ - ^™(/i, At, A;)| < m\e'^'^' - g{h, At,k)\. 

Also, 

e^'=^'-g{h,At,k) 
= e^*(^~'=-^°) - f 1 + At (jk - Jo 



= At ((4 - Jo) '{jk-Jo))+J2— {jk - JoJ , (19) 

and thus there exists Ci{h) and C2 (using ([7]) and 1^) such that 

je^^"^* - g{h, At, A;)| < At {Ci{h) + CsAt) . 

Here Ci{h) — )■ as /i — )■ 0, and C2 is bounded, see [1] for exact detail. Thus 
we get the following bound. 

Theorem 2. Accuracy: If J{x) G H^{1), r > ^ satisfies Al - A5 then 
there exist Ci{h) and C2 such that 

le^^^* - g{h, At, A;) I < At (Ci(/i) + CsAt) . 
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Applying Theorem [2] on the first part of the right-hand side of fllSp we 



get 



N_ 
2 



^ \{e'^.^'- - g-^')u^\' < At' iCih) + C^Atf \\u,\\\ (20) 



fc=-f+i 



Definition 1. fJl page 223] For integer k > 0, H^(2tx) is defined to be the 
closure of Cp{2Tc) under the inner product norm 



miH'' 



■ k 
.j=0 



(i)||2, 

L2 



For arbitrary real s > 0, H'^{2n) can also be obtained following |2l po^ges, 
219-i 



Theorem 3. JH page 223] For s G M, H^(2ti) is the set of all series 

oo 

ip{x) = ^ am^rn{x) 
m=—oo 

for which 

II II 2 I |2 I \ ^ I 1 2s I |2 ^ 

\m\*,s = l«o| + 2_^ \m\ \am\ < oo 

|m|>0 



where 



i>m{x) 



^ 



m = 0,±l,±2, 



Moreover, the norm | Iv?] |*,s is equivalent to the standard Sobolev norm \\(p\\}^s 
for^e H'{2tx). 

For exact details of the Theorem [3] please see [H page 223]. In this 
section, to use norm definitions in a periodic domain [0, 1], we use H'/'IIh^ = 

\\'^\\h'>(i)- 
Now 



^u°°^\2Tx{k + sN)) 



<^\ir^\2Tx{k + sN))\ 



< J^\u°°^^{2T:{k + sN)){2'K{k + sN)Y\^ \^\{2n{k + sN))\- 



-2a 



11 



using the Cauchy-Schwartz inequality. Thus, 



N_ 
2 

E 



N_ 
2 



^™(/i, At, k) J2 ^"'° (27r(A; + sN)) 

sjto 



9 N 

< E 









gives 



E 



^'"(/i, At, A;) J2 ^°°'° (27r(A: + sN)) 



s^O 



< /i2-C(a) Y^ \27rkf'' \u^ {2iTk)\ 

\k\>f 

< C{a)h'^\\uo\\'^C{a)h'^\\uorH.n 



(1) 



(21) 



for some a > ^. Thus applying above bounds. Lemma [3] and Theorem [3] in 
( lT6ll one gets 



|«(x,-,t^)-/7ff, < At2(Ci(/i) + C2At)'||«o||' + C3/i'"hoirH.(i)+ E 



Mfc 



|fc|>f 



< At^ {Ci{h) + C2Aty \\uor + Cah'^^lluoWi, 



■(1) 



^ ^ |fc|>f 



< At" (Ciih) + CsAt)" \\uor + Cg/i'^holll, 



■(1)' 



(22) 



for some a > ^. Thus we conclude this discussion with the following result. 

Theorem 4. // the approximation ^ of the initial value problem (J\) is 
stable, J, Jk, — Y + 1 < k < ^ satisfy Al - A5 and uq G H"{1) with (J > \, 
then there exist constants Ci{h), C2, C3(a) such that 

\\u{x,tj - f/f (x)|| < At {C^{h) + CsAt) ||^o|| + C73(a)/i^||«o||H^(i). 
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5 Convergence analysis for non-smooth ini- 
tial data 



Now let us consider an initial function which is not smooth enough so 
that II Wo II //'^(i) is bounded when i^ > |, but there exists ui < u such that 
||'"o||i^''i(i) < oo and there exists a such that 



N_ 
2 



Y^ \u°^'^{2Ti{k + sN))\ <d.x^"C{uQ). 



(23) 



fc=-T + l 



For example of such functions please see [H [15] and references there in. We 
start varying (TT6l) as 

2 



\U[X 



.t^)-Uf\\l < Y. M^y"'^ - 9'''Ul 



\k\<M 



Y kfc(0)e^'=*- - (7™[/° + Y H^)^' 



M<\k\<-^ 



o^kU, 



\k\>^ 



for some < M < y. Now with the same algebraic operations as we 
performed for the Theorem H] 



\u[x 



„U-^f||^ < 5^ |^i.(0)e^''*--r%'°| + E |^;^(0)e^'=*'"-r^^ 



39fetm _ ^m- 00,0 I 2 



|fc|<M 



M<\k\<-q- 



+ Y |^fc(0)e^'=*-|'+ Y 



|fc|>f 



|fc|<Af 



^m^^oo,0(27r(A; + siV)) 



s^Q 



E 



M<|fc|<-f 






5^ |ii,(0)e«^*- - r% 'T < t'(C^i(/^) + C2At)2|ho||', 

|fc|<M 



and 



Y i^fc(o)e^'=*'»-r^rT+ E \^^^^y 



.Qi-n 



M<\k\<i 



|fc|>- 



c E l^"'T < 



fc|>A/ 






|fc|>A/ 
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hold. Also following [H [T2] , there exists a G M such that 



< 



E 

|fc|<A'/ 

E 






|fc|<f-l 






00,0 



2n{k + sN) 
2L 



+ E 


r)^^°° 


° (27r(A; + sA^)) 


A/<|A:|<f 

9 


ST^O 




< Az2"C(wo), 





gives 



Thus we finish this section with the following theorem. 

Theorem 5. // the initial value problem (Qp with some non-smooth initial 
function Uo{x,0) = g{x) is considered so that the condition ^£3) holds and 
is approximated by the stable one step finite difference formula (EP, J, Jk, 
-f + 1 < A; < f satisfy Al - A5, then 

||M(x„U-f/f II < t(Ci(/i)+C2At)||Mo||+C2(z/i)Ax"M|Mo||f/^i(i)+Ax"C3M. 
where Ci{h), C2 and C3 are constants. 

6 Convergence of semidiscrete approximation 



Applying the discrete Fourier transform to ([ 

d ~ 

-TL^k — QkUk 

at 



where qk = [Jk — Jo) and its solution can be found as 



(24) 



Uk{t) = e^^'UkiO). 



(25) 



Applying the inverse of the DFT in 

N_ 
2 

U{x,,t)= Yl ^'''Uk{0) 



ik2TTXj 



(26) 



fc=-f+i 



14 



Thus comparing ( TT31) and (12^ 



UiXn 



oo 


t ik2nxj 


k-- 


2 

E 

JV 


1 


f/fc(0) 






N 
2 


















E 


ik2-KXj 
-1 


(uki^h 


Jkt 


- e^'''*L^fc(0)) + 


lfel>f 


,(0)e^^- 


ik27vx 



(27) 



Applying Parseval's relation to (j2 



JV 
2 



m X 



•,,t)-/7(2:,,t)|' < Yl \M^)e'''' - e^^'UkiO) + J^ |«fc(0)e'^'=*e 



Ak2iTXi 



k=-^~l 

N_ 
2 



|fc|>^ 



< Y. 1^0)6'"' - e^^'UkiO) + 5^ |«fc(0)< 

|fc|>^ 



,9*:* I 



(28) 



fe=-T-l 



Now applying Poisson's formula to the first part of the right-hand side of 



N 
2 



JV 
2 



Y \uk{0)e''' - e'^'^'UkiO) = Yl \e'''-e^''\ \u^{27rk)\' 



k=-^+i 



JV 
2 



E Y^"''''\^o{Mk + sN))[ 



fc = -f + 1 '^T^O 



Now 



Now 



Jikt _ ^qut\2 ^ |g9,t ^^ _ gffe-gfc)^ |2 < ^2|^^ _ ~^|2_ 



gfe-gfe = I Y ^k+m.N \ - Qk, using ([7]) 

v7n=— oo 
/ ^ Qk+mN = ^^ Mfc+mAf ~ -'m.iV j , siuce Qk = Jk — Jq, 

Y (</°°(27r(A; + mA^)) - J°°(27rmA)') , using (^. 
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So from [1] it follows that 

\qk-qk\<2J°°{£j=C^{h), 

and Ci(/i) — )■ as /i — )■ when J°° is smooth enough (if J°° G -^2(1^)5 then 
|J°°(OI -> as 1^1 ^ cx) [H, [ISl page 30]). Thus 

So applying the relation (TT7|) . Parseval's relation and Theorem [3] 

^ |(e^'=*-e^''=*)no°°f<t^Ci2(/i) ^ 1^^(0)1^ 

lfcl<f lfcl<4 

< t'clih) [K(o)p+5^|fcn^i.(o)f 

Thus, similar to Section H] there exists a constant 6*2(0") such that 
\\u{xj,t) - U{x,,t)f < t^Cl{h)\\uof + C^{cr)e''\\uo\\l.(^^y 
Thus we conclude 

Theorem 6. // the semidiscrete approximation ^ of the initial value prob- 
lem HP is stable J, Jk, -f + 1 < k < ^ satisfy Al - A5 and uq e H"{1) 
with 0" > |, then there exist constants Ci(/i), C2(cr) such that 

\\u{xj,t)-U{xj,t)\\ <tCi(/i)||no||+C2(a)/i"||Mob-(i). 

7 Numerical Experiments and discussion 

Here we start by experimenting numerical error in the approximation ([6]) of 
([1]). We compute numerical error at t = 1. Figure [2] shows the behaviour 
of ||it(-,t„) — U'^{-)\\h for various choices of h and At with smooth and 
non-smooth uq for all x E [0, 1] where J°°(x) = \/^exp{—10x'^). From 
Figure [21 we observe that for smooth uq the rate of convergence of the 
solutions are faster than for the non-smooth uq- Here we also notice that 
choices of h and At have an impact on the rate of convergence. That 
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Figure 2: Here we plot the error term \\u{-, tn) — U"'{-)\\h at t = 1 by varying 
At and /i with J°°(x) = -7^6"^°^' and (a) uo{x) = v/^e-^^"^)' (left figure), 
(b) Uo{x) = ie~''^'~2l (right figure). 



agrees with our theoretical estimates (and also motivates us to investigate 
further the theoretical stability and convergence analysis of such a one step 
approximation). From this study we notice that the full discrete scheme is 
conditionally stable. The accuracy of the scheme depends on the smoothness 
of the initial function. We have some limitations in this study. We impose 
some reasonable restrictions on the kernel function to prove the stability and 
convergence results. The analysis of higher order schemes in one and multi- 
dimensions leaves as a future study which is of course more challenging. 
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